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^ ■ We show that an infinite number of non-unitary minimal models may describe 

two dimensional turbulent magnetohydrodynamics (MHD), both in the presence and 
absence of the Alf'ven effect. We argue that the existence of a critical dynamical index 
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results in the Alf'ven effect or equivelently the equipartition of energy. We show that 
there are an infinite number of conserved quantities in 2D — MHD turbulent systems 
both in the limit of vanishing the viscocities and in force free case. In the force free 
case, using the non-unitary minimal model M27 we derive the correlation functions 
for the velocity stream function and magnetic flux function. Generalising this simple 
model we find the exponents of the energy spectrum in the inertial range for a class of 
conformal field theories. 
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1 - Introduction 

There has been some recent activity towards the application of Conformal Field Theory 
(CFT) to the theory of turbulence in two dimensions [1-7]. The main point is that the 
energy spectrum and higher correlation functins can be derived by means of some non- 
unitary minimal model of CFT . Polyakov [1,2] derives a few criteria for a CFT which can 
be a possible candidate for describing turbulence, and finds a candidate CFT which gives the 
value of for the exponent of the energy spectrum (experimental results give an exponent 
between 3 and 4 [11-14] ). Expanding on Polyakov's method others [3,6] argue that there 
are a large number of CFT's which satisfy Polyakov's constraints but have more than one 
primary field. The role of extra primary fields is not clear, but it has been suggested that 
they may have to do with passive scalar and magnetic fields [4]. Briefly Polyakov's method 
is as follows. To describe turbulent behaviour at high Reynold's numbers, we interpret 
the Navier-Stokes equation from a statistical mechanical point view. That is we consider 
the correlation functions of the stream function with respect to some stationary probability 
density: 

< &i 1 (x 1 )Qv t (x 2 )...&i N (x N ) > . 

as correlators of primary fields of some CFT . The requirement that these correlations be 
stationary is the main prerequisite. This condition is refered to as the Hopf equation [9]. 
The fact that the correlation functions have a power law behaviour, indicates that conformal 
invariance may be at work . Indeed in the limit of vanishing viscosity (/i — > 0) the Navier- 
Stokes equations are scale invariant ( see [8] for more details ). Also in a plasma one can 
use hydrodynamic equations to describe the collective degrees of freedom, however the more 
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conventional method uses the kinetic equations to describe turbulent processes [10]. As in 
the case of pure hydrodynamics, the two dimensional magnetohydrodynamic system (2D- 
MHD) differs from its three dimensional version in that it has an infinite number of conserved 
quantities in the limit of zero viscocity /i — > and zero molecular resistivity rj — > and in this 
limit the equations are scale invariant. Therefore 2D-MHD is also a good candidate where 
CFT may be applicable. Such a system was first considered by Ferretti and Yang [4]. Here 
we shall extend their arguments and find non- unitary models which can give the correlation 
functions and the energy spectrum index. It is worth noting that the requirement that the 
critical dynamical index be consistent, is equivalent to the Alf'ven effect, this means that 
there is equipartion of energy between kinetic and magnetic coponents and this requirement 
greatly reduces the number of possible minimal models of CFT solutions. 
This paper is organised as follows: in section 2 we first describe the equations governing 
two dimensional magnetohydrodynamic systems (2D-MHD) in the inertial range. We then 
discuss its scaling properties. In section 3 we go on to describe a conformal field theory with 
two primary fields pretaining to (2D-MHD). We then generalise to CFT's with more primary 
fields and give a table of possible solutions in the general case and in the the Alf'ven region. 
In the appendix we describe a method of deriving OPE coefficients, for any Minimal Models 
in CFT and give the coefficients up to the third level. The results are summarized in section 
4 . 



3 



2 - Properties of two-dimensional MHD 

The incompressible two dimensional magnetohydrodynamic (2D-MHD) system has two in- 
dependent dynamical variables, the velocity stream function ip, related to the velocity field 

V a 

V a = e a pdp(p (2.1) 
and the magnetic flux function tp related to the magnetic field B a via: 

B a = e a pdpip (2.2) 

here e a p is antisymmetric and e 12 = +l,e 2 i = — 1. 
The dynamics is given by the pair of equations [10] : 

— = -e a pd a ipdpw + e af3 d a ipd p J + fiV 2 w (2.3) 
dip 

-^r = -eapdaipdpip + r)J (2.4) 

where 

w = vV J = 

here \i is the molecular kinematic viscosity and r\ is the molecular resistivity. Note that 
normalisation is chosen such that the magnetic field assumes the same dimensions as velocity. 
In the inertial range rj and fi can be ignored then it follows from eqs. (2.3) and (2.4) that 
there exist three global, quadratic, conserved quantities. 

E = X - J '{V 2 + B 2 )d 2 x (2.5) 

H = J V.Bd 2 x (2.6) 



dE 
~dt 



dH 
~dt ~ 

dA 



A = J i, 2 d 2 x (2.7) 

which are the total energy, the cross helicity and the mean square magnetic potential. In 
fact this system has an infinite number of conserved quantities such as: 

R n = J *p n d 2 x (2.8) 

where n is any real number. The time evolution of the E,H and A are given as: 

-77 J {V 2 ^fd 2 x -n J (V 2 (/)) 2 d 2 x (2.9) 

-(/i + 77) /(V 2 V>)(VWx (2.10) 
dA r 

— = -r ] J (V^)Vx (2.11) 
a similar expression for the R n can be written as follows 

= _^ n ( n _ 1) J ^- 2 (V^)Vx (2.12) 

But when r\ and ji are negligible, the system of equations (2.3) and (2.4) display the scale 
invariance 

x^Xx y^A^V t^\ 1+h t V^A 1 -^ (2.13) 

Now if we impose Kolmogorov's idea of constant flux of energy we get h — | [15], so that 
the behaviour of fields under scaling are : 

x^Xx y9^A 2/3 y9 t^\ 4/ h ^-^A 2/3 ^ (2.14) 

Following Polyakov [24] we believe that this scale invariance signals conformal symmetry of 
this system. Simple scaling arguments [16-21], show that in turbulent 2D-MHD the energy 
spectrum behaves as: 

E(k)~k^ (2.15) 



Unfortunately this result is in poor agreement with recent numerical simulations [29]. But 
our results show the deviation from this spectrum and is in agreement with simulation [29]. 
Another important aspect to consider is the Alf'ven effect [22-23]. This effect essentially 
amounts to equipartition of energy between the kinetic and the magnetic components of the 
energy. Thus V fc 2 and B\ should have the same spectrum in the inertial range. We shall later 
show how this effect bears on the critical dynamical index, and limits our choice of CFT. 
If Polyakov's ideas were applicable here, a conformal field theory may exist such that its 
correlation functions coincide with those of the 2D-MHD system. However such a system 
has to be non-unitary in order to give the power law behaviour suggested by the scaling 
relations (i.e. eq.(2.14)). First of all, let us show that to describe turbulence by means of 
CFT , we have to use non-unitary minimal models. Any local conformal field Aj(x) have 
associated with it an anamolous dimension dj, i.e. under a transformation x — > \x we have 
[24,25] : 

Aj(x) -> X^Mx) (2.16) 

We therefore observe that if associated with the fields of 2D-MHD system i.e. the f and the 
if), they must have negative anamolous dimensions. On the other hand negative anamolous 
dimensions are possible only in non-unitary minimal models which in turn non-unitary min- 
imal models result in infrared divergences. We shall deal with this problem later. 

3 - Conformal 2D-MHD Turbulence. 

The simplest model we consider is M 2 j with three primary fields I , the identity </?, with 
anamolous dimension (^,^), and ip, with anamolous dimension (^, ^) . The central 
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change is C = Now we can derive the small scale behaviour of the two point functions. 

< (p( x )<p( ) >~ \x\% (3.1) 

12 

< ip( x )ip(o) >~ \x\~ (3.2) 

and the correlation function of (p and tp vanishes since they have different anamolous dimen- 
sions. However the above expressions are clearly unphysical since they grow with distance. 
To avoid this problem one has to introduce an infrared cutoff [2] . Let us look at this problem 
in the momentum space 

< tpiiQtpii-k) >~ C t \k\~ 2 - 4 ^ (3.3) 

where ipi can be either ip or tp and 

n 2 4|A W |+i r(4|A^|) 
1 r(4|A^| + 2) 

To avoid the problem of infrared divergence we can restrict ourselves to the inertial range 
i.e. i ^> K ^> , where R is the large scale boundary of the system and a is the dissipation 
range. Thus 

< (p(x)(p(0) >~ Q / k- 2 -^e lk - x dk (3.5) 

which results in 

< <p(x)<p(0) >~ R 4 ^ - x^\ + £ ( am )(-^)"Vl A ^l (3.6) 

m=l 2jK 

where 

1 r(4|A^| + l) sin(|(l-4|A^|-m)7r) 
" m m! T(2| A^| + 1 + f ) r(|(4| A^| + m + 2)) 1 ' j 

Here the IR problem just as in the case of pure turbulence [2, 3, 26], and it may be removed 
by considering of turbulent 2D-MHD with boundary. 

7 



On the other hand, there exists a dissipation scale 'a' , so that the inertial range was defined 
by A; <C K Thus fully developed turbulence is equivalent to letting 'a' tend to zero, or 
equivalently letting the Reynold's number tend to infinity. However the existance of an 
ultraviolet cutoff means that we have to be careful when products of operators at the same 
point are involved. To this end we handle such products using the point splitting technique. 
Consider the right hand side of equation (2.3,2.4) 

e a pd a ifjd(3V 2 if} = limit a ^ e al 3d a ilj{x + ^)dpV 2 i/>(x - ^) (3.8) 

e a pd a ipdf3ip = limit a ^ e a pd a (p(x + ^)dpif)(x - ^) (3.9) 

Where limit a ^ , expresses angle averaging. Now to evaluate the above, we take advantage 
of the operator product expansion. We have the general form for the fusion rule of (p, q) 
minimal model [27,28]: 

min{n\+n2 — l,2p— m— ri2 — 1) min(mi+m2 — l,2g— mi — vni— 1) 

E E (3-io) 

k=\ni— tj-2 |+1 t=\mi— ma |+1 

Where variables k, £ run over odd or even numbers if they are bounded by odd or even 
numbers respectively. For our particular choice of (q,p) = (2, 7) we have two primary fields 
<f and ip and their families [ip] and [ip] satisfy: 

MxM = [/] + M (3.11) 

MxM = M + M (3.12) 
M * [1>] = [i\ + bl>] + [v>] (3.13) 

Note that by the family of <p, we mean all operators which can be constructed from ip using 
the Virasoro generators L n , such as: 

(3.14) 
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In order to explicitly calculate the rhs. of equations (3.8) and (3.9) by means of the fusion 
rules, we use the following relations for the OPE of field operators: 

<p n (x + ^)(p m (x - |) ~ |a| 2(A ^ A "- Am) (v?p + a ia L_ lVp + a 2 (a 2 L 2 + a z L\) Vp (3.15) 

+a 3 (o 4 L 3 _ 1 + a 5 L_iL_ 2 + a 6 L^ 3 )(p p + ... 

where n,m,p = 1,2 and ipi = if, ip 2 = ip ■ In the appendix we have given a method of 
deriving the coefficients 01,02,03 etc. and have derived these coefficients up to the third 
level. By differentiation of lhs. of eq. (3.15) we will find the leading term in product of ip 
and ip in the limit a — > as follows: 

e a pd a <pd^ = Ci|a|^'(Q;L_ 1 (L_2 + sL\) - C.C)^ (3.16) 

where s is a constant, determined by the operator product expansion and o is a constant. 
It is clear that the antisymmetry of rhs. of eq. (3.16) under complex conjugation is just the 
consequence of the e-tensor at the lhs. Similar calculation shows that the other terms in 
eqs. (2.3) and (2.4) by means of OPE and point spliting procedure can be written in terms 
of Virasoro generators in the limit a — > as follows : 

e^d^dpd 2 ^ = C 2 \a\ 6 r (L_ 2 Lzl - l 2 L\)i> (3.17) 

e a pd a (pdpd 2 (p = C 4 \a\^(L_ 2 Lzl - L 2 L\)il) (3.18) 
Before proceeding further let us look at the various constants of motion R n , A 2 , if and E : 

^ = _ 1) f f-i V 2 ^ 2 I (3.19) 
at J 

For Rn, A 2 ,H and E to be constants of motion, we must have fx, rj — (77 = where /j,q and 
a are permeability and conductivity, respectively) which require the conductivity a to be 



infinite. On the other hand it is obvious that the above quantities are conserved in the limit 
of jj, and V 2 ip — > 0, that is for finite conductivity we have the same conserved quantities. 
But, letting V 2 ip tend to zero is equivalent to vanishing of J. This situation is well known 
as force free MHD [10]. Let us return to conformal field theory and look at V 2 ip — > , it is 
evident that; 

V 2 ip = 4L_ 1 Z_ 1 -0 = (3.20) 

With this condition we can proceed further with the rhs. of equations (3.16), (3.17), and 
(3.18) and we find that they are all zero. In non- unitary minimal model M2J with C = ^p, 
ip is degenarate in the second level that is: 

(L_ 2 -^L 2 _^ = 0,A^^ (3.21) 

Thus the rhs.s of eqs.(3.16),(3.17) and (3.18) vanish. The more general iV-point correlation 
function satisfies the well know BPZ equation: 

" t (I^f + < >= < 3 ' 22 > 

Now for a generelisation of our simple model let us postulate fusion rules as below: 

[<p] x M = [X x ] + ... (3.23) 

where X\ , is the primary field of the lowest dimension in the OPE of (p and ip. Equation 
(3.16) changes to 

^ ~ lim lal^-^-^iaL^ipL 2 ^ + L_ 2 ) - C.C.)X 1 (3.24) 
at o->o 

Let us now use the conserved quantities to find the constraints on dimensions of the fields. 
Similar to the arguments used in [4], we consider the cascade of mean square magnetic 
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potential: 



A 



I 



ip 2 d 2 x 



(3.25) 



this requires that < ip( x )ip( x ) >be scale invariant. 



Thus, using equations (3.24), we have 



AXi + 2 + = 



(3.26) 



AX 1 > Atp + Aip 



(3.27) 



where the second condition comes from requiring non-singularity of the rhs. of (3.24) in the 
limit a — > 0. Table 1 gives a list of models which satisfy these conditions. 
Let us now consider the implications of the Alf'ven effect [22,23]. It is well known that [30] 
the excitations of the Elsasser's fields propagate as the Alfv'en waves in opposite directions 
along the lines of force of the 'B' field at speeds of order 'B'. What is meant by the 'inertial 
range' in this context is, the region where the wavenumber of the the Alfv'en waves lie 
whitin the inertial range and the energy cascade results from the scattering of Elssaser's 
fields [30]. This implies that , there should be asymptotically exact equipartition of energy 
in the inertial range, i.e. V^ 2 = aB k 2 where a is of order unity. As discussed in detail by 
Chandrasekhar [31] , for Kolmogorov's hypotheses of similarity to holds here the constant 
a shuld take value 1.62647 . Therefore the equipartition of energy requires ip and ip to have 
similar scaling behavior: 



Aip = Aip 



(3.28) 



We can derive the dynamical index of equations (2.3) and (2.4) as : 



At = -Aip + 2 



(3.29) 
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and 

At = Ap - 2A^ + 2 (3.30) 

Where the eqs.(3.29) and (3.30) comes from eqs. (2.4) and (2.3) respectively. Thus again the 
existence of a single index for temporal scaling require Aip = Aip. We can therefore extract 
models out of the table, which are consistent with the Alf'ven effect. These are given in 
Table 2. These models are consistent with the limit // — > and V 2 t/> = 0, keeping a, finite. 
By means of the anomalous dimensions of ip and (p, we can write the energy spectrum as 
follows: 

E(k) ~ A; 4 I A ^I +1 + k 4 ^ +1 
Where the exponents of energy spectrum for are given in Tables 1 and 2. 
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4 - Concluding remarks 

In this paper we have derived a number of CFT which are possible condidates for describing 
2D-MHD. In the limit of finite conductivity requieres 2D-MHD to be force free and increases 
the number of candidate CFT's. The imposition of the Alfv'en effect is equivalent to requir- 
ing a consistent dynamical index and greatly reduces the number of candidates. A direct 
interpretation of primary fields is as yet, not possible, therefore those candidate CFT's which 
have more than two primary fields seem as plausible as M 2 j, which has two primary fields. 
Adding boundaries is the next step in this theory, which may take it closer to some physical 
problems such as the quantum wire. Work along this line is already proceeding. 
Acknowledgements. We are indebted to many of our colleagues for interesting and helpful 
discussions, in particular M.Alimohmmadi, F.Ardalan, H.Arfaei, S.Azakov and A.Yu.Morozov 
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Appendix: Calculation of OPE Coefficients 

The most general expression for the operator product expansion is [27]: 

$> n (z,z)$ m (o,o) = ^^c-M^.W^-A^^+Efe^Ap-A^-A^+Efc^W,^}^^) (A.l) 

P k 

where the coefficients are 

np,{k}~{k} _ (~<p ap,{k}np,{k} / a r>\ 

{k} = {ki, k 2 , k n } (A.3) 

Note that we have (fi(0)\0 >= \i >, for the vacuum state |0 >. Now let equation (A — 1) act 
on |0 >: 

$ n (z)|A m >= Y,C p nm z A ?- A "- A ^ p (z)\A p > (A.4) 

M^) = J2^ k ^ P ni K} L-k 1 ...L. kn (A.5) 

\z, A p >= %W > (A.6) 

Expand \z,A p > in terms of the complete basis \N, A p > where \N, A p > is defined such 
that the coefficients of expansion are z N : 

\z,A p >=J2z N \N,A p > (A.7) 

By applying of Lj over eq.(A-4) we have; 

LjlN + j^p >= (A p -A m +jA n + N)\N,A p > (A.8) 

and by solving the recursion relations we can find pP'W. For level one we have: 

Li|l,A p >= (A„- A m + A n )\A p > (A.9) 

which results is 

|1, A p >= a 1 L_ 1 |A p > (A.10) 
14 



or: 



A p A m + A n . . 

ai = 2A P (A - U) 

For the second level by means of eqs. (A-7) and (A-10) we have: 

L,|2, A p >= (A - - A " ± A " + 1)(A - - — ± K) L^A r > (A.12) 
L 2 |2,A P >= (A p - A m + 2A n )|A p > (A.13) 

which results in: 

|2,A P >= {a 2 L_ 2 + azL\)\k p > (A.14) 
where a 2 and a 3 satisfy a system of equations: 



MijCcj = At i,j = 2,3 



where: 



M = 



4A P + c/2 6Ap 



O 4Ap + 2 



A 2 




. ^ 3 . 





Ap - A m + 2A, 



(Ap-A m +A TC +1)(A P -A m +A n ) 
2Ap 



(A.15) 
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This system can now be solved to give: 



« 2 




6A P A 3 


- A 2 (AA P + 2) 


1 


«3 




3A 1 - 


A 2 (4A p + f) 


2A P (5 - 8A P ) 



A similar method will work for higher levels. For level N, we have in place of equation A-14, 
an expansion corresponding to the partition of N. We then find a system of equations by 
successively applying Lj, and finally the coefficients are derived. For the third level, using 
equation A-8, we have: 

1,13, A p >= [A p - A m + 2 + A n ]|2,A p > (A. 17) 

L 2 |3,A P >= [A p - A m + 1 + 2A„]|1,A P > (A.18) 

L 3 |3,A P >= [A p - A m + 3A n ]|A p > (A.19) 

and [3, A p > is given by: 

[3, A p >= {a A L\ + a 5 L_iL_ 2 + a 6 L_ 3 ) | A p > 

where a A ,a$ and a 6 satisfy the following system of equations: 

a 4 (24A p + 6) + « 5 {(4A P + |) + 9} + 5« 6 = B{Ap " Am) 

a 4 (24A p ) + a 5 (4(4A p + |)) + a 6 (6A p + 2c) = C 
a 4 (8(A p + 1)) + a 5 (7 + 2A P ) + 4a 6 = A(a 2 + a 3 ) 
and a 2 , a 3 are given by A-15. The coefficients A, B and C are given by: 

A = [A p - A m + 2 + A n ] (A.24) 

B = [Ap - A m + 1 + 2A n ] (A.25) 
16 



(A.20) 



(A.21) 

(A.22) 
(A.23) 



C=[A P -A m + 3A n ] (A.26) 
The inverse of the matrix of coefficients is : 

M-j 1 = (A.27) 

where 

B n = C(2A P + 3) + A P (6A P - f f ) 
B 12 = C + 3A P + 1/2 

B 13 = 1/2[C 2 + C(11A P + 8) + 2A P (12A P - 13)] 
B 21 = 4[C7(A P + 1) + 3A P (A P - 1)] 
B 22 = 2(7A P - 2) 

B 23 = 3[C(4A P + 1) + A P (12A P - 7)] 
B 31 = 2[2C(A P + 1) + 5A P (2A P - 1)] 
B 32 = 2C(A P + 1) - 8A 2 p - 38A P + 15 
£? 33 = 3[C(3A P + 1) + 2A P (12A P - 5)] 
and 

A = 2(2c 2 (A p + 1) + c(-2Aj - 16A P + 7) - A p (24Aj - 74A P + 25)) (A.28) 
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Table 1 Solution satisfying the constant mean square magnetic potential flux 
condition 
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(p,q) 


$ 




X 


4A^ + 1 


4A^ + 1 


(2,13) 


#1,2 


#1,4 


#1,5 


-2.69 


-.53 


(2, 17) 


*1,6 


#1,2 


#1,7 


-.64 


-4.88 


(2, 19) 


*1,5 


#1,2 


#1,6 


-.68 


-4.47 


(2,23) 


*1,2 


#1,3 


#1,4 


-2.30 


-.73 


(2,27) 


#1,4 


#1,2 


#1,5 


-.77 


-3.88 


(3,29) 


#2,12 


#1,3 


#2,14 


-2.17 


-1.20 


(3,34) 


#1,2 


#1,3 


#1,4 


-2.29 


-.73 


(3,37) 


#1,5 


#2,14 


#2,18 


-.18 


-5.05 


(3,40) 


#1,4 


#1,2 


#1,5 


-.77 


-3.87 


(3,46) 


#2,17 


#1,3 


#2,19 


-2.47 


-.21 


(3,50) 


#1,3 


#2,19 


#2,21 


-1.4 


-2.52 


(3,52) 


#2,21 


#1,2 


#2,22 


-.82 


-3.61 


(3,62) 


#1,5 


#2,22 


#2,26 


-.37 


-5.83 


(3, 77) 


#2,27 


#1,3 


#2,29 


-2.68 


-.36 


(3, 79) 


#1,2 


#2,29 


#2,30 


-2.10 


-.88 


(3,80) 


#2,28 


#1,3 


#2,30 


-2.7 


-.36 


(3,82) 


#1,2 


#2,30 


#2,31 


-2.10 


-.89 


(3,89) 


#2,33 


#1,2 


#2,34 


-.89 


-3.35 


(3,92) 


#2,34 


#1,2 


#2,35 


-.90 


-3.33 
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Table 2 Solutions that satisfy both the constant mean square magnetic po- 
tential and the Alf ven constraint 



(p,q) 


$ 




X 


4A„+ 1 


4A^ + 1 


At 


(4,27) 


*2,10 


*2,10 


#3,19 


-2.08 


-2.08 


1.23 


(5,32) 


*1,9 


*1,9 


#1,17 


-2.5 


-2.5 


1.12 


(5,34) 


*2,10 


*2,10 


#3,19 


-2.04 


-2.04 


1.24 


(6,35) 


*2,H 


*2,12 


#3,22 


-2.98 


-2.92 


1.0 


(6,35) 


*3,17 


*3,17 


#5,33 


-2.96 


-2.96 


1.01 


(6,37) 


*1,8 


#1,8 


#1,15 


-2.78 


-2.78 


1.06 


(6,41) 


#2,10 


#2,10 


#3,19 


-2.01 


-2.01 


1.25 
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